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ABSTRACT 



Context. Gaseous metal discs around single white dwarfs have been discovered recently. They are thought to develop from disrupted 
planetary bodies. 

Aims. Spectroscopic analyses will allow us to study the composition of extrasolar planetary material. We investigate in detail the first 
object for which a gas disc was discovered (SDSS J122859.93 + 104032.9). 

Methods. We perform non-LTE modelling of viscous gas discs by computing the detailed vertical structure and line spectra. The 
models are composed of carbon, oxygen, magnesium, silicon, calcium, and hydrogen with chemical abundances typical for Solar 
System asteroids. Line asymmetries are modelled by assuming spiral-arm and eccentric disc structures as suggested by hydrodynam- 
ical simulations. 

Results. The observed infrared Can emission triplet can be modelled with a hydrogen-deficient metal gas disc located inside of the 
tidal disruption radius, with T^g x 6 000 K and a surface mass density of E x 0.3 g/cm-. The inner radius is well constrained at about 
0.64 Rq. The line profile asymmetry can be reproduced by either a spiral-arm structure or an eccentric disc, the latter being favoured by 
its time variability behaviour Such structures, reaching from 0.64 to 1.5 R©, contain a mass of about 3-6 -10"' g, the latter equivalent 
to the mass of a 135-km diameter Solar System asteroid. 

Key words. Accretion, accretion disks - Stars: individual: SDSS J122859. 93+ 104032.9 - White dwarfs - Planetary systems 



1. Introduction 

A significant fraction of white dwarfs (WD) that have cooled 
below TefF ~ 25 000 K ( 20-30%) displays phot ospheric ab- 
sorption lines from metals (IZuckerman et al.ll2010b . These pol- 
luted WDs must actively accrete matter at rates from 10"'^ to 
lO^'^Mo/yr, otherwise the atmospheres would h ave been puri- 
fied b y gravitational settling o f heavy e lements (iPaquette et al.l 
Il986t iKoester & Wilkenll2006t lKoestejl2009h . Until the recent 
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' past it was assume d that these stars ac crete matter from the in- 
terstellar medium (iDupuis etal.111992 '). This view has changed 
considerably in the past few years because it was discovered 

\ that many of t he polluted WDs host dust discs within the stel- 
lar tidal radius (Zuckerman & Becklin'1987VBeckHn et al.'2005'; 
,Kilicetal. 2005, 2006; Farihi et al. 2007, 2009; Fai-ihi et al. 

, 12010 ). It is now widely accepted that WDs are polluted by mat- 

' ter accretion from these discs. It is thought that they contain 
material from tidally disrupted asteroids that were scattered to- 
wards the central stars as a consequence of dynamical reset- 
tling of a planetary system in the post-main sequence phase 
dDebes & Sigurdssonll2002l:lJurall2003l) . 

The photospheric metal abundance pattern in the polluted 
WDs allows us to indirectly measure the composition of the ac- 
creted matter This opens up the exciting possibility of study- 
ing the composition of extrasolar planetary material. The first 
impressive results have already been obtained, showing that 
the orbiting debris broadly mimics the terrestrial matter of 
the inner Solar Syste m, including the possibility of water 
(Zuckerman et al. 2007t iKIein et al.l 120101; iDufour et alj 1201 Ot 
[parihi et al.ii201 li) . These results are based on our knowledge of 



metal diffusion rates in WD atmospheres and envelopes, which 
are difficult to obtain (Koester 2009). Several uncertainties can 
affect the resulting composition of accreted material as con- 
cluded from the photospheric abundance pattern. One example 
is the assumption of stationary accretion on several diffusion 
timescales of all elements involved. It is therefore highly desir- 
able to exploit alternative possibilities of determining the chem- 
ical composition of the accreted material. 

Such an alternative method is offered by the recent discov- 
ery that three circumstellar dust discs around polluted WDs also 
host gaseous metal components th at are interpreted as the col- 
lision al remains of solid material (iGansicke et al.l [20061 l2007l 
l2008l) . We are developing accretion disc models to derive the 
chemical abundances in the gas discs, so that we should be able 
to directly measure the composition of the parent planetary mate- 
rial. This method is complementary to the measurement of pho- 
tospheric abundances because the disc spectra might reveal trace 
elements that are not seen in the WD spectra. It also provides 
a means to test our understanding of diffusion processes in the 
WD a tmospheres and env elopes. Building on our preliminary 
work dWemer et al.ll2009l) . we present here new results for our 
modelling efforts of gas discs. 

In Sect.|2] we briefly introduce our method, followed by a de- 
scription of the object that we study in detail (Sect.[3]l. In Sect.|4] 
we first summarise our results concerning simple models for 
pure-calcium discs, investigating the influence of effective tem- 
perature and surface-mass density on the emergent spectrum. We 
then present vertical structures and spectra of discs composed 
of an asteroid-like mixture of light metals, comprising C, O, 
Mg, Si, and Ca. Finally, we present results for non-axisymmetric 
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disc geometries in order to explain the observed asymmetry of 
the double-peaked line profiles. Based on hydrodynamical sim- 
ulations, we investigate the time evolution of this asymmetry in 
comparison with observations. 

2. Accretion-disc modelling 

For calculating geometrically thin accretion-disc models, we use 
our code AcDc (iNagel et alJl2004l) . We assume axial symmetry, 
so that we can separate the disc into concentric rings of plane- 
parallel geometry. In that way, the radiative transfer becomes a 
one-dimensional problem. By integrating the spectra of the in- 
dividual rings, we obtain a complete disc spectrum for different 
inclination angles. 

The free parameters of one ring with radius R are effective 
temperature Teff(R), surface mass density chemical com- 

position, and the WD mass Mwd- For the energy equation we 
assume that the emitted radiation is viscously generated, so the 
Reynolds number Re (or a) enters as an additional parameter In 
the case of viscous a-discs, the radial run of Tf^ffiR) and 
can be expressed in terms of the mass-accretion rate and mass 
Mwd and radius 7?wd of the central star (Shakura & Sunyaev 
1973). For comparison with observations, the emergent spectra 
from ring segments are Doppler-shifted to account for Keplerian 
motion, hence /?wd and disc inclination / appear as additional 
parameters. 

For each disc ring, the following set of coupled equations 
were solved simultaneously under the constraints of particle 
number and charge conservation: 

- radiation transfer for the specific intensity / at frequency v 
dl(v,jj,z) 



Table 1. Statistics of the model atoms used in our disc models. 
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with the absorption coefficient x, the emission coefficient rj, 
the disc height z above the midplane, and yu = cos 6, with 6 
the angle between the ray and z; 
- hydrostatic equilibrium of gravitation, gas pressure Pg^s, and 
radiation pressure 
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with p denoting the mass density and H the Eddington flux. 
Here, we also introduced the column-mass density m as 
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- energy balance between the viscously generated energy 
Emech and the radiative energy loss firad 
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with the angular velocity oj, the mean intensity 
J, and w the kinematic viscosity written following 
iLvnden- Bell&Pringlj(II974 : 



(7) 



For the models presented here we assume Re - 15 000; 
- NLTE rate equations for the population numbers n, of the 
atomic levels / 



(8) 



J*' 



where Pjj denotes the rate coefficients, consisting of radiative 
and electron collisional components. 

Detailed information about the inv olved atomic data is pro - 
vided in the form of a model atom (cf. lRauch & Deetienll2003h . 
The model atoms we used for our NLTE calculations are sum- 
marised in Table [1] They were taken from TMAD, the Tubingen 
Model Atom Databasefl 

The principal problem for any modelling attempt is posed by 
the question of what heats the Ca ii emission line region. It can- 
not be gravitational energy released through viscosity because 
the required mass-accretion rate would have to be of the order of 
10 Mo/yr, which is by many orders of magnitude larger than 
the accretion rate invoked for the p resence of settling metal s in 
DAZ photospheres lO '^ M^/vr. 'Koest er & Wilkenl l2006'). A 
speculation by Jura (2008) was additional heating by energy dis- 
sipation through disc asymmetries, which are driven b y some ex- 
ternal unseen planet. Altemativelv. lMelis et al.l (l20I0l) suggested 
a "Z n" model in analogy to H n regions. In the case of the discs, 
the metal-dominated material is photoionised and heated (hence 
the name Zii) by absorbing photons from the WD and cools 
through optically thick emission lines. Given this lack of knowl- 
edge, we need to use Tetj, which is a measure of the vertically 
integrated dissipated energy (Eq.|5), hence the amount of energy 
radiated away from the disc surface per unit time and area, as a 
free parameter 



http://astro.uni-tuebingen.de/~TMAD/TMAD.html 
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3. SDSSJ1 22859.93+1 04032.9 

Our models are tailored to SDSS J122859.93 + 104032.9 (hence- 
forth SDSSJ1228 + 1040). This metal-polluted WD was the 
first one discovered t o be surrounded by a gaseous metal disc 
(iGansicke et al.l2006l) . It is a DAZ white dwarf with atmospheric 
parameters T^b = 22 020 ± 200 K and log ^ = 8 .24 + 0.04, and the 
derived stellar mass and radius are Mwd - 0.17 + 0.02 and 
/?WD = 0.0111 +0.0003 Rq. 

The Can infrared triplet {AA 8498, 8542, 8662 A) with 
double-peak emission line profiles is the hallmark of the 
gaseous metal discs (Fig.[B. In the case of SDSS J1228+1040, 
[Oansicke et al. (2006) measured a peak-to-peak separa- 
tion of 630km/s, i.e. the Keplerian rotation velocity is 
V sin; = 315km/s. From a spectral analysis with a kinematical 
LTE emission model, it was concluded that we see a geomet- 
rically thin, optically thick gaseous disc at high inclination 
(i - 70°). Two other weak emission featur es of Pen AA 5018, 
5169 A were seen bv iGansicke et al.l (l2006l) . Subsequent obser- 
vations bv'Mehs et al] (1201 Oh failed to detect the Fe ii ^ 5018 A 
line, possibly because of the lower signal-to-noise ratio of 
their high-res olution spectra as compared to the low-resolution 
spectrum of iGansicke et al.l (12006) . On the other hand, very 
w eak emissions from the Ca ii H and K lines were discovered 
bylMelis et al.' ("2010"). 

There is a clear asymmetry in the emission strengths of the 
double-peak Une profiles of SDSS J1228+1040 dGansicke et alJ 
I2OO6I) . A similar phenomenon is well known from Be star 
discs (jCarciofi 2010) and is ascribed to one-armed spiral waves. 
In addition, Meli s et al.l (12010) observe that the asymmetry in 
SDSS J 1 228 H- 1040 had changed such that the stronger of the two 
emission peaks has switched from t he red side of the do uble- 
peaked emission complex, as seen in IGansicke et al.l (12006). to 
the blue side. We describe these characteristics in more detail in 
Sect. l4.3l where we investigate the temporal variability predicted 
by our models. 



4. Results 



4.1. Parameter study for pure calcium discs 



In a first exploratory study (IWemer et al.ll2009l ), we calculated 
disc models composed only of calcium, with two values for Z 
and three values for T^ff. The inner and outer disc radii were 
set to 1 .0 Rq and 1 .2 R©, respectively. We found that the emis- 
sion strength of the Can triplet decreases with T^ff increas- 
ing from 5 000K to 7000K (Fig.[B- The reason is the shift- 
ing Can/Cam ionisation balance. A closer comparison of the 
three line components shows that their emission strengths be- 
come equal with increasing Tetj, a behaviour that constrains Teff . 
A similar trend is seen when Z is reduced from 0.6g/cm^ to 
0.3g/cm^ at reff = 6 500K. We stress that the models have a 
considerable continuum flux compared to the line-emission peak 
heights. The relative strength of the profile depression between 
the double-peaks increases with increasing inclination, the dou- 
ble peak structures become broader The disc models are opti- 
cally thin in terms of the Rosseland optical depth. 

The observed spectrum does not show H(j, emission. This can 
be used to determine an upper limit for the hydrogen abundance. 
We varied the H content (H= 1%, 0.1%, 0.01%, by mass) and 
found that with an abundance of 1%, the H„ peak height is com- 
parable to that of the Ca 11 triplet, so would be detectable in the 
spectrum of SDSS J1228-H1040. 
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Fig. 1. Normalised spectra of three pure Ca disc models (top) 
with different Teff compared to the observed spectrum of 
SDSS1228-H1040 taken from SDSS (bottom). The model spectra 
are shifted vertically for clarity. 




log m [g cm"' 



Fig. 2. Vertical structure of every other of the 1 1 disc rings along 
the column-mass density (increasing from outer layers towards 
the disc's midplane). 



The effective temperature of the disc is well constrained by 
three models, with reffai5 800K (Fig.[T]), but the asymmetry 
of the line profiles is of course not matched by our symmet- 
ric models. The cooler model (Teff = 5 000K) is perhaps more 
favourable because of the larger line-to-continuum emission ra- 
tio, while the hotter model (Teff^VOOOK) has the advantage 
that the relative strengths of the three line components are re- 
produced better 
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Fig. 3. Spectra of three (non-rotating) circumstellar rings at 1-11° and: (a) 2.0, (b) 71.4, and (c) 136.4Rwd- For clarity, the spectra 
are vertically shifted by amounts as colour-coded numerically on the plots. With the exception of the infrared Can triplet and the 
adjacent C ii doublet, as well as Ca ii H & K, no fine structure splitting was taken into account. 



4.2. 0-Si-Mg-C-Ca-H disc models 

In the next step, we expanded the set of considered chemical ele- 
ments in order to achieve a composition comparable to CI chon- 
drites in the Solar System or a bulk-Earth mixture (iKlein et al.l 
^010). New species included are O, Si, Mg, and C. Iron poses 
special numerical problems and will be introduced in future 
work. The initially chosen element abundances are representa- 
tive of the class of CI chondrites. In detail, they are: H = 10"^, 
C = 4.6, O = 65.5, Mg = 13.5, Si = 15.1, Ca = 1.3 (% mass 
fraction). 

We investigated the influence of the radial disc extent on 
the spectrum. At its largest, the disc model consists of 11 
rings extending from R, =2Rwd = 0.022 R© to an outer radius 
of Ro - 136RwD = 1.5 Rq. The ring radii and effective tempera- 
tures are listed in Table|2] The model for ring 8 did not converge; 
its spectrum was set equal to that of ring 9. The entire disc is as- 
sumed to have a radially constant surface density of 0.3 g/cm^. 

In Fig.|2] we present the vertical run of temperature, 
Rosseland optical depth Tross, mass density, gravity, and geomet- 
rical height of the disc rings. Only the inner rings are optically 
thick. Figure|3]shows spectra of three different rings, at 2.0, 7 1 .4, 
and 136.4RwD- They all show emission lines, getting stronger in 
the outer parts of the disc. 

We combined the 11 rings to several axisymmetric discs 
with varying radial extent. Their spectra are shown in Fig.H 
The dip between the triplet line components of the Can lines 



Table 2. Radial position and effective temperature of the 1 1 con- 
centric rings forming the basic disc. 
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seen in observed spectra can best be reproduced if the disc does 
not reach too far inwards, resulting in a minimum inner radius 
of about 60RwD- Variation in the outer disc radius has almost 
no effect on the Ca ii triplet, but reducing the outer radius de- 
creases the line strength of a Cn doublet {AA 8685, 8699 A) 
near the red component of the Ca n triplet. While the lack of ob- 
served C n emission may suggest a constraint on the outer disc 
radius, we note that decreasing the carbon abundance by an or- 
der of magnitude also removes this feature from the model spec- 
tra (see Fig.|5]l. Given the carbon-poor nature of the observed 
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Fig. 4. Effect of variation in the disc's radial extent at / =77°. 
Upper panel: the inner radius is fixed at 58 Rwd and the outer ra- 
dius varies: 84, 110, and 136 Rwd- The Cii line disappears with 
decreasing outer radius. Lower panel: The outer radius is fixed 
at 1 10 Rwd and the inner radius varies: 2, 30, and 58 Rwd- Broad 
line wings develop when the inner radius becomes smaller 
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Fig. 5. Variation in the carbon abundance. Shown are two models 
(a: 4.6% and b: 0.46% by mass) and the observed spectrum (red). 
The radial disc extent is 58-136 Rwd, ' = 77°. 



discs (iReach et al.ll2005t lJura et al.ll2009l) and polluted WD at- 
mospheres ('J urall2006l) . abundance may play a substantial role. 

Our model strongly overpredicts the strength of the Ca ii 
H & K emission lines. This proble m has already been noted 
in the pure-Ca disc described above (IWerner et al.l 12009*) . This 
problem did not disappear in the present multi-element disc 
model. In contrast, our model features emission lines from other 
metals, markedly from Mgii that are not observed. Another 
weakness of our model is the relatively strong continuum flux 
that would detectable by distorting the WD spectrum. This prob- 
lem is alleviated when the eff'ective temperature is reduced. 

4.3. Asymmetric disc models 

To investigate the asymmetry of the line profiles in the spectra 
of SDSS J1228-I-1040, we modified our method combining the 
disc rings to receive a complete disc spectrum. In the surface 
integrating step, we used only parts of the rings in order to con- 
struct a spiral-arm like or an eccentric shape of the disc. The ring 
segments are still assumed to undergo Keplerian rotation for the 
calculation of the spectral Doppler shift. Depending on the ori- 
entation towards the observer, these non-axisymmetric accretion 
discs result in asymmetric line profiles that can be compared to 
the observations. 
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Fig. 6. Two exemplary structures from the hydrodynamical sim- 
ulations. The initial gas blob first evolves into a spiral-arm struc- 
ture (left) and after some time into an eccentric disc (right). The 
surface density is colour-coded as indicated. 
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Fig. 7. Sketch of the two non-axisymmetric geometries used for 
our spectral models. Left: a spiral-arm like structure. Right: an 
eccentric accretion disc. Both are composed of disc ring seg- 
ments (red). For the unmarked ring segments, the emergent flux 
is set to zero. The observer's position is towards the right. 



We performed hydrodynamic calculations with the FARGO 
code (Masset 2000) in order to motivate our particular choice 
of geometries. This code was originally developed to compute 
the hydrodynamic evolution of a protoplanetary disc on a fixed 
polar coordinate system, but it is also suitable for other kinds of 
sheared fluid discs. A possible scenario could be that an asteroid 
coming within the tidal radius of the WD is disrupted, forming 
a locally concentrated debris cloud. Matter then spirals inwards. 
We start the simulations by putting a gas blob at the tidal radius 
(R = 137 Rwd) onto a circular orbit (PqA = 3.6 • 10"^ s) around 
the WD. The gas mass is 7-10^' g. The initial surface density 
distribution is Gaussian with a blob radius of 2.8- 10^ km. We 
chose an open boundary condition and fixed the simulation rim 
at a value of R - 194 Rwd- 

The Keplerian rotating material gets smeared out into a 
spiral-arm like structure within a short time of f = 1.21-10'*s. 
After another Af = 3.1-10"'s, the spiral is catching up with its 
own starting point and an eccentric closed disc forms and re- 
tains this shape for the rest of the simulation. For both sit- 
uations the surface mass density is rather homogeneous with 
2* w 0.3 g/cm^. Examples for both geometries found by the sim- 
ulations are shown in Fig.|6] In Fig.|2]we display the correspond- 
ing assembly of our ring segment models. 

The upper panel of Fig. [8] shows the resulting spectrum in 
the case of a spiral-arm like shape with a minimum radius of 
58 Rwd and a maximum radius of 1 36 Rwd for an observer's po- 
sition according to Fig.|7](left). The strength of the Ca ii lines and 
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Fig. 8. Emergent spectrum for the (a) spiral-arm like and (b) ec- 
centric disc structures as shown in Fig.|7] The modelled profiles 
(thick lines) display an asymmetry as observed (thin line). For 
the spiral arm the radial extent is 58-136Rwd, whereas in the 
case of the eccentric disc it is 45-136 Rwd- For both, / =77°, 
= 285°, and i; = 0.3 g/cml 
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Fig. 9. Time series of the computed Ca ii triplet emission, cover- 
ing 100 min. For the spiral-arm geometry (left panel), a change 
in the line-profile asymmetry evolves very fast, whereas for the 
eccentric disc (right) it takes much longer. 

the asymmetry of the blue and red parts of each line component 
are well reproduced. Only the dip between the Doppler-shifted 
parts of each line are not as deep as in the observation. The mass 
of such a spiral arm would be about 3.9 ■ 10^' g. In the lower 
panel of Fig.[8j the resulting spectrum of an eccentric disc for an 
observer's position according to Fig. Q (right) is shown. The fit 
quality to the observed spectrum is similar to the spiral-arm case. 
Such an eccentric disc would have a mass of about 6.3 • 10^' g. 

4.4. Time variability 

iGansicke et al.l (l2006h present two time series of spectra from 
SDSS J1228H-1040 taken on June 30 and July 1, 2006. They in- 
dicate that the asymmetry of the Ca ii line profiles did not change 
during the two 100 min observations. To investigate this with our 
models, we assumed a Kepleiian rotating spiral arm and an ec- 
centric disc with (jj = 0.01°/s, and calculated synthetic spectra in 
steps of = 6° concerning the orientation towards the observer, 
starting at = 231° and centred on 285°, for which we found the 
best fit to the emission lines. The resulting spectral time series 
are shown in Fig.|9] covering 100 min i n order to be compara - 
ble with the lower right part of Fig. 1 in IGansicke et al.l (l2006h . 
In the case of a spiral-arm like geometry, the asymmetry of the 
line profile would change significantly within this time intei-val 
(Fig.m left panel). The separation of the maxima of the double 
peaks decreases, and at the same time the Can line asymme- 



try reverses in wavelength. In contrast, for the eccentric disc the 
relative strength of the asymmetry changes quite slowly (Fig.|9] 
right) and the double-peak separation remains almost constant, 
which is in better ag reement with the n on-variable observations. 
On the other hand, iMelis et al] (12010 ) found a switch of the 
asymmetry on a longer timescale, between the IGansicke et al.l 
observations in 2006 and their own in 2008. 



5. Summary and conclusion 

We performed non-LTE modelling of gaseous metal discs 
around WDs in order to study their spectral signatures in com- 
parison to observations of SDSS J1228-I-1040. The modelling 
was done in three steps. 

At first, pure calcium models were constracted to constrain 
the disc characteristics by fitting the observed infrared Ca ii 
emission triplet. Qualitatively good fits can be obtained with a 
geometrically and optically thin, Keplerian viscous gas disc ring 
at a distance of 1.2 from the WD, with Teff ~ 5 800 K and a 
low surface mass density 2"si0.3g/cm^. The disc is hydrogen- 
deficient (H < 1% by mass), and it is located within the tidal dis- 
ruption radius (/?tidai = 1.5 R©). 

In the second step, we constructed axisymmetric disc models 
composed of elements in Chondiitic abundance, namely C, O, 
Si, Mg, and Ca. We found that the inner radius of the observed, 
emitting, Ca ii gaseous component of the disc is well constrained 
at > 0.65 Rq = 58 Rwd- The outer radius can be constrained by the 
emission strength of a C ii doublet {AA 8685, 8699 A) that is not 
seen in the observations. An alternative explanation could be a 
reduced carbon abundance, which would be a hint that the disc 
in SDSS J1228+1040 has a bulk-Earth like composition instead 
of a CI chondritic one. 

In the third step, we investigated asymmetric disc structures 
by assuming spiral-arm and eccentric disc shapes as suggested 
by hydrodynamical simulations. Both geometries can qualita- 
tively explain the asymmetry observed in the double-peak line 
profiles well. An investigation of the time variability of the 
computed line profiles suggests that the eccentric disc model 
displays less significant variability than the spiral-arm geome- 
try. Considering the current observational material, the eccentric 
disc model is more realistic. Mass estimates for the circumstel- 
lar gas material using the two geometric models results in 3- 
6- 10^1 g, which are equivalent to the mass of a 135-km diameter 
Solar System asteroid. 

Acknowledgements. We thank Tobias Miiller from the Computational Physics 
group at our institute for providing us with the FARGO simulations. S.H. and 
T.R. are supported by DFG (grant Wel312/37-l) and DLR (grant 50OR0806), 
respectively. 

Note added in proof. Following the acceptance of our paper, 
an HST/COS observation of SDSS J1228 + 1040 became pubHc 
(dataset LB5Z040, observation date 2010-02-19). This is the first 
available UV observation. The spectrum shows a strong emis- 
sion line of Mg ii that our disk model, as displayed in Fig. [3] pre- 
dicts to be the strongest UV line. This discovery is essential be- 
cause Mg is the third element (after Ca and Fe) that is observed 
in the metal disk, supporting the idea of ground planetary mate- 
rial. Figure [To] shows the observed Mg ii resonance doublet { AA 
2796.35, 2803.53 A) in emission with central absorption compo- 
nents of photospheric origin. Overplotted is a pure-H WD model 
scaled to fit the continuum and our disk model (including all 
rings 5-11 as listed in Table |2] / = 77°) scaled arbitrarily, plus 
the co-added WD-i-disk model spectrum. The optical depth in the 
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Fig. 10. Mg II resonance doublet in the HST/COS spectrum, 
overplotted with our (rotating) disk model and a pure-H WD 
model, plus the co-added model. Fine-structure splitting of the 
line is considered in the disk model. 

Mg II line varies across the disk and has a maximum of r ~ 10^ 
in the outermost disk region. 
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